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Comment on "Density of States and Critical Be- 
havior of the Coulomb Glass" 

In a recent Letter pj, Surer e£ a/, concluded that 
their simulation results are consistent with the Efros 
Shklovskii prediction for the density of states in the three- 
dimensional Coulomb glass. Here, we show that this 
statement has no relevance concerning the problem of the 
asymptotic behavior in the Coulomb gap since it is based 
on unjustified assumptions. Moreover, for the random- 
displacement Coulomb glass model, we demonstrate that 
a part of the density of states data by Surer et al. er- 
roneously exhibit a broad gap. This is related to the 
staggered occupation being instable contrary to [l| . 

In detail, Efros and Shklovskii considered systems of 
localized charges (see the review |2|). They showed that 
the single-particle density of states p vanishes when the 
energy E approaches the chemical potential, which equals 
here. For stability with respect to arbitrary single- 
particle hops, they analytically derived that, as E — > 0, 
p(E) = a\E\ s with S = 2 in the three-dimensional case. 

Surer et al. [l] state correctly, that the literature sit- 
uation concerning quantitative numerical studies of this 
Coulomb gap is not yet satisfactory. Then they claim 
that their p(E) data "can be fit very well with a form 
oc \E\ S " with 5 close to 2. This is puzzling for their sam- 
ples are considerably smaller than the samples in our old 
work [3j], where clear deviations from oc E 2 were found. 

Inspection of Fig. 1 of [1] uncovers a first problem: The 
simulation data seem to be approximated by the ansatz 
p(E) — a\E\ A +b instead of by the Efros Shklovskii power 
law, but this is not mentioned in the text of 

For a detailed analysis, we digitized Fig. 1 of [l[ by 
means of WinDig. The data for the disorder strength 
W = 0.4 are re-plotted in a log-log representation in 
Fig. 1(a) here. Their precision is confirmed by a fit to 
p{E) = a\E\ s + b for \E\ < 0.3 and L = 14 yielding 
6 = 1.84(3), in almost perfect agreement with [l[. 

Figure 1(a) shows that these p(E) data clearly do not 
follow a pure power law, so that the quality of the fit 
could be ensured only by the unphysical constant b. 

One might argue, b emulates the finite-size effects not 
analyzed in [l[ , being important due to the long range of 
the interaction Q. However, the data of [l[ exhibit an 
unexpected feature contradicting this idea: For sample- 
edge length L = 14, b is clearly larger than for L = 8, 
although b should vanish as L — ► oo [4j. A second con 
arises from the p(E) by Surer et al. being far larger than 
the data from Fig. 3(b) in [3j for small \E\, albeit another 
random potential was used there (compare 0]). 

To clarify these points, we performed own simulations 
for W = 0.4 and various L using periodic boundary con- 
ditions and minimum image convention. The results are 
included in Fig. 1(a). For large \E\, our p(E) obtained by 
relaxation concerning all single-particle hops nicely agree 
with p(E) from [lj, but for \E\ < 0.3, within the whole fit 
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FIG. 1: (Color online) p(E) for Gaussian random-potential 
(a) and random-displacement (b) Coulomb glass models with 
W — 0.4 and 0.2, respectively, (a) From [if: y, L = 8, and 
A, L = 14. Own: Y, L = 6, x, L = 20, and +, L = 70, 
for complete relaxation via single-particle hops, whereas □, 
L — 6, and 0, L = 20, for relaxation only via short-range hops 
(see text), (b) From [J: continuous line. Own: Q, charged- 
ordered occupation without relaxation, +, relaxed concerning 
single-particle hops. All three for L — 14. Error bars of own 
data are omitted since they do not exceed the symbol size. 


region of [l[, they are considerably steeper and smaller. 

According to the original reasoning of the Coulomb gap 
2] , this difference may arise from the minimum search in 
l| not considering all single-particle hops. Figure 1(a) 
includes two p(E) curves which we obtained by relaxation 
only via short hops ranging up to third-next neighbors 
(length < V3). These relations support our hypothesis. 

The situation is even more inconsistent in case of the 
random-displacement model: Compared to [l[, we found 
far stronger broadening of the 6 peaks of a charge-ordered 
(NaCl) occupation of the lattice. Figure 1(b) shows that 
already for W = 0.2, a weak Gaussian random displace- 
ment, broadening causes p(Q) » 0.11. Hence, a Coulomb 
gap is formed by relaxation, in which the perfect charge 
order claimed in Fig. 2 of [l( for T — > is disturbed. 
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